Dephasing by time-dependent random potentials 
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Diffusion of electrons in a two-dimensional system with time-dependent random potentials is in- 
vestigated numerically. The correction to the conductivity due to inelastic scatterings by oscillating 
potentials is shown to be a universal function of the frequency u>, which is consistent with the weak 

2 

localization prediction % logo;. 



Transport properties of the two-dimensional disordered 
electron systems have attracted much attention in re- 
cent years. To understand the transport properties of 
a random system, the concept of quantum interference 
plays an important role. The constructive interference 
of two time-reversed trajectories enhances the probabil- 
ity for backscattering, leading to the contribution called 
weak localization correction. Rj3| Another example is the 
phenomenon of the sample specific reproducible conduc- 
tance fluctuation, p,pj which is highly sensitive to the 
motion of a single impurity atom in the metallic region. 
p-pl More recently, such impurity-motion-induced con- 
ductance fluctuations have been utilized experimentally 
to study the detailed dynamics of impurity atom tunnel- 
ing in a single two-level system in disordered metals. 0] 
The motion of an interstitial impurity between two en- 
ergetically equivalent stability points that have different 
effects on the conductivity is also studied to understand 
the temporal fluctuations of the conductance, where one 
can observe 1// noise. 10 

Such effects of quantum interference are suppressed in 
the presence of the inelastic scattering, for instance, due 
to the electron-phonon interaction or due to the Coulomb 
interaction of conduction electrons. There exists a time 
scale called dephasing time t^, after which the interfer- 
ence effects are destroyed. fhylJJt The dephasing time 
is determined by the inelastic process and usually de- 
pends on the magnetic field or on the temperature. Con- 
sequently, the effect of interference can be observed ex- 
perimentally by controlling these parameters . 

In this paper, we study the effect of dynamically fluc- 
tuating potentials on the conductivity by solving the 
time-dependent Schrodinger equation numerically. We 
adopt here the formula for solving the time-dependent 
Schrodinger equation proposed by Suzuki, p3J-fi8[ which 
make it possible to deal numerically with the time- 
dependent potential on the sufficiently large system of 
two-dimensional lattice. We observe the second moment 
of the wave packet, which gives information about exten- 
sion of the wave packet. We then evaluate the conduc- 
tivity and discuss how the conductivity depends on the 



motion of potentials, in particular, on the frequency of 
potentials. 

In order to solve numerically the time-dependent 
Schrodinger equations, we adopt the method based on 
the higher-order decomposition of exponential operators. 
|l3| , [l9| The basic formula we have used is the fourth-order 
decomposition of exponential operators: 



(1) 



■■e xAl/2 , 

(2) 

and the parameter p is given by p — (2 — \/2)~ l . Here 
A\, . . . , A q are arbitrary operators. 

We consider the tight-binding Hamiltonian with time- 
dependent potential £i(t) on the two-dimensional square 
lattice: 





exp[x(Ai + A 2 + ■■■ + A q )} 




= S(xp)S (x(l - 2p)) S(xp) + 0(x 5 ) 


where 




S(x) = 


e xA 1 /2 e xA 2 /2 _ _ _ xA q _i_/2 xA q xA g _ 



H(t) = 



(y) 



iWJCi 



(3) 



where Cj(Ci) denotes a creation (annihilation) operator 
of an electron at the site i and (i,j) is the nearest neigh- 
bor site. We then decompose this Hamiltonian into five 



parts as described in ref. 



H(t) = Y,Hn 



n=l 



H X = V 



H 2 = V 



E 
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E 
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namely, 



H 5 (t), 



z2 c f+x c p- 



(4) 



h.c. 
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Cl + -C P + h.c. 



H A 



vj: 



e 
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c-t a-+h.c. 



H 5 (t) = Y,^{t)ClC P , 



(5) 



where x(y) denotes the unit vector in the x(y)— direction 
and r x {r v ) the x(y) component of the position vector r 
of the sites. All the length-scales are measured in units 
of the lattice constant clq . Each Hamiltonian H n consists 
of operators which commute with each other, hence we 
can obtain the analytical expressions for exp(—iH n t/K) 
by diagonalizing two by two matrices. 

The state vector \ip(t + St)) at time t + St is obtained 
as 

m + st)) = u(t + st,t)\m), (6) 

where the time-development operator U(t + St, t) is de- 
fined by 



U(t + St, t) = T exp 



t + st 



H(s)ds 



(7) 



with T the time-ordering operator. Using the formula 
proposed by Suzuki, Jlj this ordered exponential can be 
expressed by an ordinary exponential operator as 



U(t + St,t) = exp 



St(--H(t) 



T) 



where the super-operator T is defined by 



F(t)e dU G(t) = F(t + St)G(t). 



(8) 



(9) 



Using formula ([!]) and (||), the time-development opera- 
tor (pi) is decomposed as the product of exponential op- 
erators; p7| 

U(t + St,t) = S 2 (-iStp/h, t + (l-p/2)St) 
x S 2 (-i5t(l - 2p)/h, t + St/2) 
x S 2 (~iStp/h, t+pSt/2) 
+ 0(St 5 ) (10) 



with 



S 2 (x,t) 



2 xH 1 / 2q xH 2 /2 _ , e xHi/ 2e xH 5 (t) 



^xHa/2 . 



, e xH 2 /2 e xH!/ 2 



(11) 



In order to consider the time evolution of the wave 
packet with fixed energy E, we have carried out a nu- 
merical diagonalization of a subsystem Hn whose size 
is No by Nq, located at the center of the whole system 
and have chosen an eigenstate of -ff/v with eigenvalue 
En sa E as the initial wave packet. 

The quantity we observe is the second moment of the 
wave packet (r 2 (t)) c defined by 



(r 2 (t)) c ^(\r(t)\ 2 )-(mY' 



(12) 



2 W> 



dr2r d - 1 drr 2 |V'(f,i)|' 



(13) 



where ^(r,t) denotes the wave function at time t, and d 
the dimensionality of the system. If the wave function ex- 
tends throughout the whole system, the second moment 
is expected to grow in proportion to time t 



= 2dDt. 



(14) 



Here, the coefficient D denotes the diffusion coefficient. 
In contrast, if the wave function is localized, it is clear 
that the second moment remains finite in the limit t — > 
oo. Pjl In the metallic region the Einstein relation 



a = e 2 p(E F )D 



(15) 



relates the conductivity a to the diffusion constant, where 
p(E-p) is the density of state at Fermi energy E-p. 

In order to take into account the effect of moving po- 
tential, we assume that the site-potentials take the form: 



£? {t) = e?(0)cos[ojt + e(r)}, 



(16) 



where u> is the frequency. Effects of scattering from impu- 
rities are introduced through randomness of site energy 
ep(Q) and phase 6(f) at t = distributed uniformly in 
the regions, 



MO) | < W/2, 
\6(r)\<n. 



(17) 



We consider the adiabatic case u> <ti V/fi, where impuri- 
ties move slower than electrons. 

We have calculated the second moment of the wave 
packet (r 2 ) c at various random potential strength. The 
size of the systems are 500 by 500 for the energies 
E-p/V = — 1. We have carried out an exact diag- 
onalization for the 20 by 20 subsystem at the center 
of the system and taken the eigenfunction of the sub- 
system whose eigenvalue is closest to the given energy 
E-p as the initial wave packet. By this procedure we 
can simulate the diffusion of the wave packet whose en- 
ergy is approximately equal to Ep. The single time 
step St is taken to be Q.2(h/V) in the simulation. With 
this choice of St, fluctuations of the expectation value 
of the Hamiltonian (H) = (ip(t)\H\ip(t)) for ui = can 
be safely neglected [19100] throughout our simulations 
(t < 2000(ft/V)). 
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FIG. 1. 

Fig. [H gives an example of the calculated second 
moment for uj = lCT 3 V/?i, W/V = iy/2. The 
second moment is proportional to the time in a wide 
range up to 2 x 10 3 h/V. We can evaluate the diffu- 
sion constant of this system from Eq. (fL4j) with least 
square fit to these data and obtain the diffusion con- 
stant of the sample to be D = lAQa 2 V/Ti. In the 
actual simulation, the quantities (r 2 ) c are averaged over 
at least five samples of random potential distribution. 
The density of states p(E-p) in Eq. (Il5|) is evaluated 
from the direct diagonalization of the two-dimensional 
lattice model for 40 by 40 sites for u> — 0, and an av- 
erage over 100 samples is performed. We have obtained 
that p(E F )alV = 0.292, 0.291, 0.281, and 0.251 for 
W/V = 2, 2\/2, 4, and 4a/2, respectively. 
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In Fig. ||, we show the conductance for lo — ob- 
tained by the present method, and compare it with that 
calculated by a conventional recursive Green's function 
method. [£lj The results show that two methods agree 
fairly well with each other. We can also see the con- 
ductivity is proportional to (W/V)~ 2 in a wide range of 
W/V. 
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In Fig. [| we show examples of the conductivity as the 
function of lnw for W/V = 2^2,4, and 4\/2, where 
lu = Lu/(V/h). The conductivity linearly depends on 



In a) as 



AIiiuj 



(7(0). 



(18) 



The coefficients A(« 0.22) calculated by the least square 
fit to the averaged conductivity, are almost the same 
among these three different disorder cases, although the 
conductivity <r(0) for u> = is significantly different 
from each other. 

To interpret the above universal correction due to fre- 
quency w, we recall the weak localization correction, 



5a = 



2n 2 h 



■ ma. 



(19) 



where r^ is the dephasing time, and r the elastic scatter- 
ing time. The lo dependence of t^ in the kicked rotator 
problem |22| as well as that in the one-dimensional dis- 
ordered system j23[ in the presence of noise is estimated 



to be 

5a M 



u> 



' 2 / 3 . Then the weak localization correction 



is estimated to be 



5a 



— -In — 

37T h LOq 



(20) 



The pre-factor 2/3tt — 0.212 ... of the lnw term agrees 
with the numerical calculation shown in Fig. 0, which is 
universal and independent of the potential strength W. 
Note that this ^-dependence is also observed in the case 
of low energy phonon scattering. |2q| 

To justify our argument, the condition r <§C r^ is nec- 
essary. This is why the conductivity deviates from the 
logo; behavior in large w region, where r is the order of 
t^. We also see the deviation from the logw behavior in a 
small u> region. In this case, w is too small for the poten- 
tial to be changed in a finite time where our simulation 
has been performed. 

In conclusion, we have analyzed the weak localiza- 
tion effect in the two-dimensional system with time- 
dependent random potentials using the equation-of- 
motion method. It has been shown numerically that the 
weak localization correction to the conductivity due to 
the fluctuating potentials takes universal value indepen- 
dent of the random potential strength W. The correction 

9 2 

term agrees with the formula 5a = -^ ^- In ^- , proposed 
for the inelastic scattering by moving impurities. Our re- 
sults will open a new way to incorporate the dephasing 
mechanism in numerical simulations. 
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Fig. |l|: An example of the second moment as the function 
of time for w = lQ- 3 V/h,W/V = 4^/2, and E F /V = 
— 1. The solid line corresponds to D = lA6a 2 V/h. 
Fig. |2j: Conductivity obtained by the equation-of-motion 
method (solid line) for uj = 0. The dashed line corre- 
sponds to the conductance calculated by the recursive 
Green's function method. 

Fig. g: The calculated conductivity of the two- 
dimensional system with time dependent impurity po- 
tential for (i) W/V = 2^2, (ii)W/V = 4, and 
(iii) W/V = 4V2. The solid lines indicate the cor- 
rections da 
5a = 0.214 



0.232^- In w, 5a 



0.217^- lnw, and 



■Law, respectively. 



